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Abstract 

We address the controversial hot question concerning the validity of the loose coupling versus the 
lever-arm theories in the actomyosin dynamics by re-interpreting and extending the phenomenological 
washboard potential model proposed by some of us in a previous paper. In this new model a Brownian 
motion harnessing thermal energy is assumed to co-exist with the deterministic swing of the lever- 
arm, to yield an excellent fit of the set of data obtained by some of us on the sliding of Myosin 
II heads on immobilized actin filaments under various load conditions. Our theoretical arguments 
are complemented by accurate numerical simulations, and the robustness of the model is tested via 
different choices of parameters and potential profiles. 
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1 Introduction 

In 1985 an experimental set-up to measure the distance traveled by an actin filament interacting with 
a Myosin II filament 1 during a complete ATP cycle was described [1]. Under low load conditions on 
the actin filament, traveled distances of up to 60 nm were observed. This value was not understandable 
according to the most popular idea of the sliding mechanism, rotation or tilting of the myosin head 
bound on actin molecule, conceived as directly coupled with the ATP hydrolysis (see [2]). Indipcndent, 
successive experiments (see [3], [4] and references therein) then showed distances ranging between 4 and 
10 nanometers. Such smaller range of traveled distances is coherent with the lever-arm theory (see [5] 
and [6]) that is still widely believed to account for the generation of the force responsible for the actin 

*Work performed within a joint cooperation agreement between Japan Science and Technology Corporation (JST) and 
Universita di Napoli Federico II, under partial support by MIUR (PRIN 2003) and by Campania Region. 
1 Hereafter we shall use the word "myosin" as an abbreviation of "Myosin II" . 
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filament sliding. The level-arm theory can be essentially schematized as follows (more details can be 
found in [7], where a graphical animation is also present): 

- The binding of an ATP molecule on the catalytic site of myosin head produces the detachment of 
the head from the actin filament. 

- The ATP hydrolysis, that follows within about the next 10 ms, changes ATP to ADP.Pi, determining 
a free energy variation in one reaction cycle of about 20 fc B T(see [8]). ATP hydrolysis also generates 
a sensible rotation of the neck of the myosin head. 

- The transition M. ATP — > M. ADP. Pi is accompanied by an increase of affinity of this complex 
for actin (see [7]), which enhances the chance of myosin to get binded on an actin site. Since such 
binding is reversible, the myosin head can visit more than a single actin site. 

- While myosin head is binded at an actin site, its neck may suddenly switch back to its original 
position (this is the so called "power stroke" ) with the successive release of the phosphoric radical 
Pi and the occurrence of conformational changes around the myosin coiled coil 2 . This process is 
reversible i.e. phosphoric radical Pi can re-establish the complex M. ADP. Pi. In the absence of 
such recombination, the above-mentioned conformational changes are very likely to generate the 
sliding of the actin filament with the release of ADP molecule after about 2 ms. 

- The myosin head stays then binded on the actin site till the arrival of a new ATP molecule, which 
starts the whole process afresh. 

Summing up, as far as reference is made exclusively to the mechanical phenomenon of the generation 
of the force responsible for the movement, the lever-arm theory is strictly deterministic: Each ATP cycle 
generates one single power-stroke that causes a sliding of constant length and preassigned direction of 
the actin filament. Within such framework, a tight coupling between ATP cycles and protein movements 
is envisaged. 

Even though the apparent disagreement of the above mentioned difference of travelled distances could 
be understood on the base of a low duty-ratio characterizing Myosin II heads (see [10], pages 221-226), 
successive experiments by Kitamura et al. ([11]) raised anew the question of such a disagreement. Indeed, 
aiming at a continuous efforts towards increasingly accurate measurements of the traveled distances during 
an ATP cycle, it was possible to achieve great improvements of the technological set-up that included 
design and construction of "home made" highly sophisticated devices. By exploiting such a technology, 
it was proceeded as follows: A single myosin head was attached to the tip of a glass microneedle and 
placed near an actin filament that had been previously immobilized on a microscopy slide by means of 
optical tweezers. The deflections of the needle with respect to its resting position were then measured and 
recorded. From the obtained traces the following three features emerged that are in evident disagreement 
with the tight coupling theory: 

1. The total traveled distance (i.e. the total displacement) of the myosin head is not constant and it 
can be as large as 30 nm. 

2. This displacement is the sum of a random number of single "steps", the amplitude of each of which 
equals the distance (~ 5.3 nm) between two successive actin monomers. During the time elapsing 
between two successive steps, the myosin head randomly jitters around an equilibrium position. 

3. Steps mainly occurs in a fixed "forward" direction, although some of them occasionally take place in 
the opposite "backward" direction. Hereafter, forward steps will be taken as positive and backward 
steps as negative. The total displacement is thus the algebraic sum of the number of performed 
forward and backward steps. 

Such evidence contrasts the one-to-one relation hypothesized between ATP hydrolysis and the occurrence 
of the mechanical event consisting of the power stroke in the myosin head. Furthermore the observation 
of the existence of random elements leads one to conjecture that a significant role could be played in 

2 The nature of such conformational changes has not been yet fully understood. With reference to kinesin, a possible 
mechanism is described in [9] where its validity also for the myosin family is conjectured. 
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such context by the thermal agitation of the environmental molecules of the watery solution in which the 
involved proteins are embedded. This is the motivation for the assumption of the existence of a loose 
coupling between ATP cycle and actomyosin dynamics (See, for instance, [2]). 

While referring to [12] for a lucid outline of the origin of the controversy existing among the supporters 
of the loose coupling approach and the community of those faithful to the lever-arm theory, and therefore 
also to the tight coupling vision, in the present paper we purpose to extend the phcnomenological model 
earlier proposed by some of us ([13]) by including in it the effect of the swing of the lever-arm. Thus 
doing we show that the experimental data of myosin sliding obtained in [1 L] under various load conditions 
can be very well accounted for. 

Specifically, via the comparison of our theoretical results with those yielded by our experiments, we 
shall test the assumption that during the time interval elapsing between the attachment of myosin head to 
the actin filament and the final release of the phosphoric radical (rising phase in the sequel), the position 
of the myosin head is determined, both by the lever-arm swing and by the action of Brownian motion, 
that includes a macroscopic deterministic force responsible for a non-zero average net displacement. Such 
a direction-orienting force can be viewed as the originated by changes of chemical states that accur within 
the myosin head and that are fuelled by energy supplied by myosin head itself. Such a view is coherent 
with the notion of "effective driving potential" in the sense, for instance, of Wang and Oster [14]. As a 
consequence, the coupling between the ATP cycle and the mechanical effects should appear to be less 
rigid, and the actomyosin dynamics should definitely exhibit variability features of the kind mentioned 
in the above 1.-^3. items. 

The above outlined conjectured random dynamics will be described in details in Section 2 with 
a specific reference to an earlier paper [13]. Here we limit ourselves to showing that, within such a 
framework, it is possible to handle a key point of the mentioned controversy. Indeed, in [12] the essential 
relevance of the role played by the length of the myosin neck is stressed, since the lever-arm theory makes 
the sliding distance less than, and somewhat proportional to, such length. This would for instance imply 
that reducing the length of the myosin neck to a half should reduce to a half the sliding distance, which 
appears to be confirmed by the experiments in [15]. On the contrary, certain performed experiments only 
show slight changes of the overall displacement of the myosin head even after complete removal of its 
neck [12]. 

In order to reconcile these evident discrepancies, we shall assume that the displacement X of the 
myosin head during the ATP cycle could in general be envisaged as a linear combination as follows: 

X = rX R + dX D (1) 

where Xr denotes the displacement induced by the biased thermal effects and Xjj the displacement 
generated by the power stroke, and where r and d are constants, each of which hereafter will be taken 
as equal to or 1. Note that the case r = and d = 1 yields the lever-arm theory; instead, setting 
r = 1 and d = depicts a purely random situation, in the absence of any sliding due to the power stroke. 
Finally, the case r = d = 1 leads to an integration of the two theories. With the choice r = d = 1, 
the controversial results related to the length of the neck of myosin head can be overcome by assuming 
that the random displacement is a few times larger than the deterministic one. This assumption implies 
that only slight changes of the total distance traveled by myosin head would be observable when the 
contribution Xo due the power stroke is made somewhat smaller by shortening the length of the myosin 
neck. 

2 The model 

Let L denote the distance between each pair of neighboring actin monomers. As suggested in recent 
literature (see [1G] and [17]) in our computations we shall take L = 5.5 nm. In addition, we shall 
assumed that the magnitude of the sliding induced by the swing of the lever-arm equals that of a step, 
and thus take Xjj ~ L. 

Our conjectured relation (1) with r = d = 1 is preliminarily supported by the bar chart in Fig. 4c) 
of [11] showing that at least one step is performed by each and every myosin head during the rising phase. 
The second column of Table 1 shows the heights of the columns of the mentioned bar chart, whereas the 
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first column indicates net step numbers, i.e. the integral part of the ratios X/L. The quoted experiment 
was performed under low load conditions by means of microneedles having stiffness less than 0.1 pN/nm. 

Table 1: Observed distribution of net number of steps performed by myosin 
heads as indicated in [11]. Here [X/L\ is the total step number during the 
entire rising phase. 



[X/L\ 


Observed frequency 


[X R /L\ 


Theoretical frequency 


1 


14 





15 


2 


21 


1 


22 


3 


18 


2 


17 


4 


10 


3 


8 


5 


3 


4 


3 


6 





5 


1 


total 


66 


total 


66 



As shown by the first two columns of Table 1, one sees for instance that out of 66 observed myosin 
heads (all attached to a glass microneedle of stiffness less than 0.1 pN/nm) 14 performed a unit net 
steps, 21 a net step number equal to 2, etc., to conclude that 3 of them have performed a net step 
number equal to 5 implying total displacement of about 30 nm. Columns 3 and 4 of Table 1 show that 
[Xr/L\ = \_X/L\ — 1 is well fitted by a Poisson distribution with parameter fj = 1.5 given by the ratio 
of the total number of performed net steps (99) to the number of considered myosin heads (66). 

The agreement between experimentally observed frequencies and those predicted via the Poisson dis- 
tribution, jointly with the rarity of the backward steps, leads one to conclude that the "dwell time" 
(namely the time interval elapsing between two successive steps) should be, to a good approximation, ex- 
ponentially distributed. Such conclusion is experimentally supported by the data summarized in Fig. 4b) 
of [11] leading to the estimate of about 5 ms for the mean dwell time. 

Our facing sequences of rare events with exponentially distributed intcrarrival times is strongly sug- 
gestive of a first-exit problem out of an interval for continuous Markov processes possessing an equilibrium 
point sufficiently far from at least one of the end points of the diffusion interval (See, for instance, [18]). 

On the ground of all foregoing considerations, with reference to the random part of the rising phase 3 
we are led to construct a model for the actomyosin dynamics that is based on the following assumptions: 

(i) The complex M.ADP.Pi + Energy is viewed as a point-size particle moving along an axis X on 
which the abscissa x denotes the displacement of the particle from the starting position. The 
positive direction of X is that of the forward steps of myosin head. 

(ii) The particle is embedded in a fluid. Hence, it is not only subject to a dissipative viscous force 
characterized by a drag coefficient /?, but also to microscopic forces originating from the thermal 
motion of the fluid molecules. On account of the fluctuation-dissipation theorem, such microscopic 
forces can be macroscopically described by means of a Gaussian white noise having intensity 2{3k B T, 
where k B is Boltzmann constant and T the absolute temperature. 

(iii) The global interaction of the particle with the actin filament is synthesized in a conservative unique 
force deriving from a potential U(x). The structure of the actin filament suggests that U(x) be a 
periodic function with period equal to the distance L between pairs of consecutive actin monomers: 

U(x) = U(x + rL), VreZ. (2) 

Henceforth we shall denote by La (0 < La < L) the minimum of U(x), assume U(La) = and 
denote by Uq := U(0) = U(L) the depth of the potential well, namely the maximum of U(x). 

3 By "random part of the rising phase"we denote the time interval during which the myosin head is subject exclusively 
to Brownian motion. 
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(iv) The particle's dynamics is described by Newton's equation in which the total acting force is the 
sum of two terms: The first term is a deterministic force generated by the potential U{x) and by 
the viscous force, while the second term is the random force due to the presence of the Gaussian 
white noise. 

(v) Two more (constant) forces, denoted by Fj and F e , act on the particle. We assume that F{ is 
generated by a process that finds its origin in a part of the energy possessed by the particle, and 
take FiL <C Uq. Instead, F e is an external force, conceivable applied from the outside by the 
experimenter. 

Summing up, we are assuming that the complex M.ADP.Pi + Energy can be looked at as a Brownian 
particle subject to a tilted potential V(x): 

V{x) = U{x) - Fx (3) 

where U(x) has been indicated in (2) and F = Fi — F e . In the experimental conditions the height of 
the potential wells is Uq < 100 pN-nm, the period of the potential is L — 5.5 nm, the particle mass 
is m = 2.2 ■ 10~ 22 kg, the drag coefficient is (3 = 90 pN-ns/nm and the environmental temperature is 
T = 293 K. Therefore, Reynolds's number is much less than 1 (see, also, [19]), so that the inertial term 
of the equation of motion can be disregarded. In conclusion, the overdamped equation describing the 
movement of the particle is the following Langevin equation: 

2KT A(t) (4) 

where A(t) is a zero-mean white Gaussian noise with unit intensity, ( ' ) denotes time derivative and ( ' ) 
space derivative. 

Within this framework, this idealized Brownian particle randomly moves around an equilibrium point 
located at a minimum La of U(x). Whenever it exits the current potential well, we conventionally 
say that the corresponding myosin head has made a step, in the forward or in the backward direction 
according to where the exit has taken place. Hence, we arc facing a first-exit problem of the Brownian 
particle from the endpoint of the current potential well. Taking into account the above assumptions, we 
conclude that the distribution of the first-exit time of the Brownian particle from the potential well is 
exponential, since 

- the process described by equations (2), (3) and (4) possesses an equilibrium point at the minimum 
La of U(x); 




- the time for the particle to travel the distance La in the presence of the only force due to U (x) is 
(see for instance, [20]) r = [3L 2 a /Uq\ 

- the standard deviation of the Gaussian steady-state distribution of the process modeling the parti- 
cle's motion is a = y/ (k B T/ (3) ■ r/2 = La/ \/2uq~, where we have set uq = Uo/k B T; 

- the ratio I a '■= La/<t = \/2uq falls well within the interval (2,4-\/2) for all choices of uq S [2, 16] to 
which we shall refer in the foregoing, so that the exponential approximation for the first-exit time 
is valid [18]. 

In the sequel, we shall exploit the known formulas [21] for the conditional probability p that the 
particle moves a step forward after leaving the current potential well, 

P = l + cxp{-FL/k B T)' (5) 
and for the mean first-exit time \i from a potential well: 



;i = dxcxp feU-/ yex n-^T (6) 
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3 Determination of Fi and Uq 

Hereafter, we shall view as constants the parameters f3 and T that characterize the thermal bath and the 
period L of the potential U(x), their values being given in Section 2. Hence, the quantitative specification 
of the model described by Eqs. (2), (3) and (4) requires that numerical values be attributed to three more 
parameters: the depth Uq of the potential well, the position La of the minimum of U(x) in (0, L) and the 
internal force Fi. The numerical specification of these parameters can be performed after the function 
U(x) has been chosen. We shall preliminarily take U(x) as a symmetric (La = L/2) saw-tooth potential 
(see Table 6). Successively, we shall test the robustness of our model by assuming alternative potential 
functions, henceforth called "potential profiles". It should be noted that Fi is somewhat related to the 
largest force that myosin is able to endogenously generate. Here, we shall not attempt to provide any 
biological justification of its genesis. We limit ourselves to pointing out that elsewhere [22], where similar 
experiments were performed on Myosin VI, it is conjectured that the week binding between actin and 
myosin is a source of distortion of the geometry of the two helixes in the actin filament. Such a distortion 
exposes the hydrophobic region of actin to myosin head, thus generating a tilt of the potential, and hence 
a constant force Fi. The existence of a tilt of the potential, conjectured in [13], has been successively 
supported by means of simulations in [23] where it is shown that in the absence of such a tilt experimental 
available evidence on the myosin motion in the presence of contrasting applied loads cannot be accounted 
for by any of the other models therein considered. Within our strictly phenomenological framework the 
existence of this force Fi is supported by Eq. (5) showing that in the absence of external applied forces 
(i.e. F e = 0) steps on either directions would be equally likely unless Fi ^ 0, in contrast with the 
experimental evidence on the high degree of directionality exhibited by motion of the myosin head. 

To proceed along the quantitative specification of our model in a way to be able to attempt the fitting 
of available experimental data, the values of Fi and of Uq must be specified. This will be done by making 
use of the available experimental data [24] shown in Table 2 and in Table 3 4 . From them it is evident that 



Table 2: For three conditions of the applied load 
C, measured at the end of the rising phases, the table lists 
the recorded numbers n/ of forward steps, fib of backward 
steps, the total number of steps and the percentage p of 
forward steps. 



Table 3: Recorded dwell times fi 
for different values of the applied 
load C. 



C (pN) 


% 


fib 


hf + m 


P= 

rif + rib 


] 0.0, 0.5] 


54 


9 


63 


0.8571 


] 0.5, 1.0] 


40 


9 


49 


0.8163 


] 1.0, 2.0] 


29 


19 


48 


0.6042 



C ( P N) 


fi (ms) 


0.046 


5.3 


0.190 


5.7 


0.300 


6.0 


0.470 


7.1 


0.690 


8.9 


0.830 


6.2 


1.240 


11.1 


1.890 


11.0 



the frequency p of forward steps decreases as the applied load increases, while the dwell times increase 
with the load. 

For some fixed values of internal force Fi, use of Equation (5) has been made to calculate the theoretical 
probabilities of the particle's exit from the current potential well to the next well (i.e. the analogue of 
the forward step probabilities) as function of F e . The results are shown in Figure 1, where eight realistic 
values of Fi have been chosen in the interval [l.OO pN, 1.90 pN] . Vertical lines indicate the 3 load 
intervals of Table 2, whereas horizontal lines indicate the 3 corresponding recorded frequencies. We see 

4 Note that loads and dwell times in Tables 2 and 3 must be viewed as averages to which confidence intervals are 
associated. For instance, the load 0.046 pN is the result of all measurements for which the product of the stiffness of the 
glass microneedle times the distance traveled by the myosin head falls around 0.046 pN. The corresponding dwell time 
5.3 ms is to be viewed as the arithmetic average of the dwell times recorded during these measurements. 
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that for internal forces Fi = 1.00 pN and Fi = 1.90 pN the plotted curves do not meet the requirement of 
leading to the experimentally recorded frequency p = 0.8571, whatever value F e is chosen within interval 
[0 pN, 0.50 pN]. Similarly, for Fi = 1.55 pN no value of the computed probability equals the frequency 
p = 0.8163 for F e ranging in [0.50 pN, 1.00 pN] . Instead, all remaining 5 curves referring to values of Fi 
ranging from 1.60 pN to 1.80 pN in steps of magnitude 0.05 pN. from below upward, are in agreement 
with the experimental values of Table 2. Hence, the interval of values for Fi to be selected in order to 
secure the fitting of the experimental data is [l.60 pN, 1.80 pN] . 

1.0 -i 1 




0.3 - 



0.2 -I 1 1 1 ^ 

0.0 0.5 1.0 1.5 2.0 

External applied force (pN) 

Figure 1 : The conditional probability p that the particle moves a step forward after 
leaving the current potential well is plotted as function of the external applied force 
F e for various values of internal force Fi (pN) . Temperature and potential period are 
taken as T = 293 K and L = 5.5 nm, respectively. The three horizontal lines whose 
ordinates are 0.8571, 0.8163 and 0.6042, respectively, denote the sample percentages 
p of forward steps under the three load conditions indicated in Table 2. 

We now come to the estimation of the last parameter, Uq, i.e. of the depth of the potential well. This 
is done by exploiting Eq. (6) in which the left-hand side is viewed as the function fj, = fJ,(Fi, F e , Uq; L, /?, T) 
and fixed to 5.3 ms, that (see Table 3) corresponds to the smallest applied load 0.046 pN. Since temper- 
ature T. period L of the potential V(x) = U(x) — (Fi — F e )x, and drag coefficient (3 arc specified, if we 
take F e = 0.046 pN Eq. (6) makes Uq an implicit function of Fi. 

Note that as Fi increases (i.e. as the potential tilt increases) the mean first-exit time \i decreases. 
Hence, in order to keep /x constantly equal to 5.3 ns, while Fi ranges in the interval [l.60 pN, 1.80 pN] , 
depth Uq must be taken as a monotonically increasing function of Fi. From /x(1.60, 0.046, Uq) = 5.3 and 
/i(L80, 0.046, Uo) = 5.3, Eq. (6) yields Uq « 15.632 k B T and U w 15.755 fc B T, respectively. 

In order to test the agreement of our model with the experimental values of mean dwell times for 
the various loads (see Table 3), we make use of Eq. (6) to determine /x as a function of F e for the pairs 
(1.60 pN,15.632 k B T) and (1.80 pN,15.755 k B T), involving the extrema of the determined values for Fi 
and Uq- The obtained values are showing Figure 2, where the corresponding experimental mean dwell 
times are also indicated. Obviously, all other pairs of admissible values of Fi and Uq lead to curves lying 
inside of the above two pairs. Inspection of Figure 2, jointly with the magnitudes of the confidence 
intervals associated to each experimental value [17], shows that for small values of F e the agreement 
between experimental data and theoretical predictions is good. Indeed, up to the first 3 values of Table 2 
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12.0 -r 

11.0 - 



5* 10.0 - 




5.0 H 1 1 1 1 1 H 1 1 1 — 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

External applied force (pN) 



Figure 2: First-exit time fi is plotted as a function of the external applied force 
F e . The chosen values of Fi (pN) and Uo (k B T) are indicated for each curve. Dots 
represent the experimental dwell times of Table 2. Drag coefficient (3 has been taken 
as 90 pNns/nm, whereas temperature and period of the potential are the same as in 
Figure 1. 

the agreement is excellent, to become more than acceptable for larger values of F e up to 0.83 pN. The 
discrepancies shown by the remaining two experimental data are a consequence of the crude assumption 
(vi) above by which the applied load is assumed to be strictly parallel to the direction of motion. This is 
acceptable for reasonably small loads, but certainly unrealistic for large loads: In this case a significant 
orthogonal component should be expected, whose effect is alike to an increase of the depth of the potential 
well Uq. Nevertheless, in the interval between the last two load values of Table 2, a qualitative similar 
behavior of theoretical curves and experimental values is present. 

4 Net step number distribution 

In the present Section we shall implement the model described by Eqs.(2)-(4) in order to obtain, via a 
suitable simulation procedure, the distribution of the net number of steps performed by the myosin head 
during the time interval elapsing between the instant when the ATP molecule is hydrolyzed and the final 
release of the phosphoric radical, i.e. during the random part of the rising phase. The results of our 
simulations will then be discussed with reference to the experimental data of Table 1. We recall that 
such data refer to 66 myosin heads that have globally performed 99 net steps. 

Our simulation procedure is based on a discretized version of Eq. (4) and on a routine for generating 
Gaussian pseudorandom numbers in a way to determine step by step the positions achieved by 66 Brow- 
nian particles all originating at xq at time 0. The traveled distances at the end of the individual random 
rising phases are recorded, which finally leads one to the determination of the 66 net step numbers. Such 
procedure has been implemented on an IBM SP4 parallel supercomputer and repeated 1600 times to 
obtain a reliable statistics. 

In order to solve Eq. (4) numerically, one has to specify preliminarily initial position xq of the Brownian 
particle and duration of each sample path. We shall safely take xq = L/2 since, whatever the actual 
initial position of the particle inside of the potential well, the its relaxation time is much smaller than the 
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mean first-exit time. Incidentally, we note that such a choice is consistent with the procedure adopted in 
Kitamura et al. [11] (see legend to Figure 2d), where to determine the average recorded trajectory, the 
starting positions of the rising phases have been synchronized. The specification of sample path duration 
is somewhat more involved because is a random variable of unknown distribution. Nevertheless, we can 
appeal to the approximate exponential distribution of the first-exit time, motivated by the depth of the 
potential well, to assume that © is approximately gamma-distributed, with probability density r(/i, v) 
where /i is the mean first-exit time from a potential well and v is the mean total exits of the Brownian 
particle. While [i is obtained via Eq. (6), our estimate v of v is obtained by using the data of Table 1 
and Eq. (5): 

- - 99 

V ~ 66(2p- 1)' 

Finally, for each specified potential well Uq, the size At of the time parsing in Eq. (6) is determined 
by progressively reducing it until the obtained distribution becomes appreciably invariant. For an im- 
mediate comparison, Table 4 shows the distribution obtained via simulation for a potential Uq — 5 k B T 
and a parsing time At = 0.25 ns jointly with the experimental distribution. The agreement between 
observed and numerical frequencies is more than satisfactory. The unique somewhat apparently appre- 
ciable discrepancy is that 2 (out of 66) simulated Brownian particles are seen to perform 1 negative net 
step number, implying a zero total displacement due to the superposition of the effect of the final power 
stroke. However, at the present stage, such a discrepancy should not be taken as particularly significant. 
Indeed, the presently available experimental setup does not allow one to monitor simultaneously the tra- 
jectory of the myosin head and the occurrence of hydrolysis of an ATP molecule. Therefore, observing at 
least one positive net step number of myosin head indicates that one ATP hydrolysis has occurred, and 
consequently such a case has been included among the recorded data. Instead, in the case when a zero 
net step number is observed, it is impossible to claim that an ATP molecule has been hydrolyzcd, and 
hence such cases have been disregarded. Future endeavours could aim at the design of an experimental 
setup able to prove that a zero net step number can actually occur as a consequence of the hydrolysis of 
an ATP molecule. 

Table 4: Net step frequency distributions for Uo = 5 k B T and for three values of Fi chosen 
in the admissible interval. Here F e = 0, xo = L/2 and At = 0.05 ns. Other parameters have 
been chosen as in Figure 2. The last 5 rows show the coefficient of variation, probability of 
a step forward and mean exit time, all obtained via simulations. The theoretical values of p 
and /x, obtained via Eqs. (5) and (6), are also indicated. 



[Xr/L\ 


Observed frequency 


Numerical frequency 






Fi = 1.70 pN 


F, = 1.75 pN 


Fi = 1.80 pN 


-1 




2.05 


1.90 


1.76 





14 


13.79 


13.63 


13.88 


1 


21 


20.72 


20.83 


20.93 


2 


18 


16.41 


16.43 


16.47 


3 


10 


8.53 


8.59 


8.54 


4 


3 


3.14 


3.26 


3.16 


5 





0.93 


0.93 


0.88 


CV 


0.966 


0.964 


0.963 


P 


0.909 


0.914 


0.920 


P 


0.910 


0.915 


0.920 


-p (ns) 


1235.057 


1207.288 


1180.852 


H (ns) 


1233.784 


1206.048 


1178.695 



The chosen value of Uq is motivated by a twofold consideration. First of all, it is large enough to 
imply that the first-exit time distribution is approximately exponential. This is indeed supported by the 
numerically evaluated variation coefficient that has been found to be 0.96. Second, and most relevant 
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consideration, is that it is reasonable to conceive that, under the assumed overdamped regime, the net 
step number is insensitive to the depth of the potential well. Indeed, the forward exit probability p 
given by (5), under the fixed environmental temperature, only depends on the energy FL, namely on the 
difference of potential V(x) over one period, thus being independent of the depth Uq. Table 5 evidently 
supports such conclusion. Indeed, it indicates that, for instance, by doubling the depth of the potential 
well Uq, the net step frequencies are not affected significantly, even for different choices of internal forces 5 . 



Table 5: Net step frequency distribution numerically obtained. Here F e = and xq = L/2. Other 
parameters are chosen as in Figure 2. Parsing steps are chosen as follows: At — 0.1 ns for Uq = 4 k B T, 
and At = 0.025 ns for Uq = 8 k B T. The indicated values of Fi belong to the admissible interval. 





Fi = 1.70 pN 


Fi = 1.75 pN 


Fi = 1.80 pN 


[X R /L\ 


Uq = 4 k B T 


U = 8 k B T 


Uq = 4 k B T 


U = 8 k B T 


U = 4 k B T 


Uq = 8 k B T 


-1 


2.11 


2.05 


1.94 


1.90 


1.85 


1.81 





13.56 


13.96 


13.54 


13.76 


13.48 


14.00 


1 


21.13 


20.26 


20.92 


20.26 


21.30 


20.46 


2 


16.70 


16.11 


16.93 


16.29 


16.98 


16.09 


3 


8.43 


8.58 


8.53 


8.61 


8.40 


8.56 


4 


2.95 


3.40 


3.01 


3.53 


2.91 


3.49 


5 


0.75 


1.12 


0.78 


1.15 


0.76 


1.11 



Note that the parsing parameter At must be determined with specific reference to the magnitude of Uq. 
Indeed, At must be much smaller than the relaxation time of the force —U'{x) oc j3L 2 /Uq. Furthermore, 
it must be such as to cope with the random forces due to the thermal bath. Therefore, the mean square 
displacement per unit time should remain constant as the depth Uq of U (x) is made to change, which 
implies an inverse dependence of the magnitude of At on the square of magnitude of the potential well. 
This is the motivation for the choices of At indicated in Table 5. 

5 The role of potential forms and asymmetries 

In [13] the idea and role of a washboard-potential played in myosin dynamics was introduced and exploited 
with the only reference to a parabolic potential (parabolic "profile" , in the currently adopted terminology) , 
though without any consideration on model robustness. This task is now accomplished in the more general 
framework of the present model, by considering not only parabolic, but also other types of profiles: saw- 
toot, cosine- like and Lindner- type (see Table 6). The effects of the natural asymmetry present in the 
acto-myosin system is investigated within our model with asymmetric saw-toot potential. A plot of the 
considered potentials over one period are shown in Figure 3, whereas Figure 4 refers to the corresponding 
generated forces. Lindner-type potential has been indicated for two values of parameter 8. It is not 
difficult to see that 8 — ► yields the cosine profile, whereas the potential flattens down in the middle as 
(5 increases. 

The independence of the exit probability p of the potential's profile implies that the net step distri- 
bution is profile-independent as well. This is also evident from Table 8, where the net step distributions 
are reported for each of the above-considered four potential profiles. 

Next task is to pinpoint the effects of the potential profiles on the mean first-exit time. To this 
purpose, we refer to Table 3 showing that the mean dwell time in lowest load condition is fj, = 5.3 ms. 
After choosing Fi = 1.75 pN, we then make use of Eq.(6) for each and every one of the four considered 
potentials imposing that the left hand side equals such value of /u. By iterated numerical integrations, 

5 The presence of non-integer numbers in Tables 4, 5 and 8 is a consequence of our adopted estimation procedure. The 
half-width of the related 95%-confidence intervals has been seen never to exceed 0.2. For comparison reasons we have not 
rounded out the raw numbers to the nearest integers. 
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Tabic 6: Potentials' profiles. 



Saw-toot 



U s (x) 



-Up 
U 



L 



(x — La), < x < La 
— (x-L A ), L A <x<L 



Parabolic 



Up(x) 



£ 2 /4 



L 

x 

2 



Cosine 



Uc{x) = ^ 



2tt 



Lindr 



U L (x) 



Up 

=,2.5 _ 



5 



i 27r 
cos I —X ) + 1 



- 1 



Table 7: For each potential pro- 
file the depth Uq of the poten- 
tial well is indicated. Here Fi = 
1.75 pN, F e = 0.046 pN and the 
other parameters are the same as 
in Figure 2. 



Tvnp 

j i 


Tin (k T) 

U P V^B^ ) 


Saw-toot 


15.723 


Parabolic 


15.043 


Cosine 


13.944 


Lindner (5 = 0.1) 


13.942 


Lindner (5 = 0.5) 


13.918 


Lindner ((5 = 1) 


13.851 


Lindner (S = 2) 


13.697 


Lindner (5 = 3) 


13.611 


Lindner (5 = 4) 


13.591 


Lindner (5 = 5) 


13.604 


Lindner (5 = 10) 


13.749 



for each potential profile the corresponding value of Uo is finally determined. The result are listed in 
Table 7, where the effect of the parameter 5 in Lindner-type profile has been detailed. 

The behavior of mean first-exit time /j, as a function of the magnitude of the external force is shown 
in Figure 5. All considered cases of Table 7 lead to graphs falling in the region bounded by the lowest 
and the highest curves. We point out that changing the potential profiles never yields mean first-exit 
time changes exceeding 10%. Hence, we are led to conclude that the depth Uo of the potential well can 
be tuned to acceptable biological values by a suitable selection of the potential profile, without affecting 
appreciably the value of the mean first-exit time. For instance, switching from saw-toot to Lindner-type 
potential with 5 = 2, lowers Uo by more than 2 k B T. (See Table 7). It is thus conceivable that a 
variety of potential profiles exists such that the height of the potential well can be further lowered, in 
a way to switch from about 15 k B T as indicated in [23] to about 5 k B T as suggested in [20]. The 
quantitative analysis and comparison with available data has been performed under the assumption of 
rigorous symmetries exhibited by the profiles of the potentials generating the periodic force acting on 
the Brownian particles. It is, however, presumable that the complex biological reality underlying the 
observed motion of the myosin heads may require to relax the rigorous symmetry assumption. To test 
the effect of symmetry breaking, we take into consideration the saw-toot potential Us(x) of Table 6 and 
make it asymmetric by taking La ^ L/2, namely by shifting in either direction the point of minimum of 
the potential. Thus doing, the introduced asymmetry should affect the motion of the Brownian particles. 
While the probabilities of exit from the current potential well are insensitive to the potential's profile, 
and hence also to its asymmetry, the mean first-exit time [i is clearly affected by it, as shown by Eq. (6). 
The quantitative dependence of fj, on the potential's asymmetry is indicated in Figure 6, where on the 
abscissa the external force F e is indicated. Each curve is characterized by two parameters: the point La 
of minimum and the depth Uo of the potential. From top to bottom, the first curve is characterized by 
asymmetry La = 2.0 nm, implying a somewhat steeper rise of the potential in the backward direction. 
Such asymmetry is mitigated in the next curve (La = 2.5 nm). The third curve, indicated in bold, is 
used as a comparison tool, as it refer to symmetric profile (La = L/2). The remaining two curves are 
labeled by La = 3.0 nm and La = 3.5 nm, thus representing mirror-image situation with respect to the 
zero-asymmetry case. Now the forward edges of the potentials are steeper. For each curve, the indicated 
value of Uo has been determined by imposing that /i = 5.3 ms when F e = 0.046 pN (lowest load), so that 
all curves originate at the same point. As Figure 6 shows, changes of the potential's asymmetry of 30% 
in either direction do not affect greatly the magnitude of the mean first-exit time. Finally, we remark 
that suitable choices of backward asymmetry (such La = 2.0 nm in Figure 6) improve the fitting of the 
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Figure 3: Plots of the potentials as function 
of the position within each pair of consecutive 
monomers. 
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Figure 4: Conservative forces originated by the 
potentials of Figure 3. 



experimental data even for larger load values. 



6 Concluding remarks 

The model of actomyosin dynamics discussed in the foregoing rests on the assumption that the total 
energy made available to the myosin head by the ATP molecule hydrolysis and by the thermal bath has a 
two-fold overall role: To produce the power stroke predicted by the lever-arm model and also to generate 
the kind of sliding of myosin head on the actin filament in the "loose coupling" mechanism originally 
hypothesized in [25] and then experimentally demonstrated in [11]. This is expressed mathematically via 
the representation of the displacement of a particle consisting of a combination of a deterministic part 
and of a random component. The latter is generated by the simultaneous presence of a washboard- type 
potential and a random force arising from thermal fluctuations. Our model has then been tested by 
making use of a set of data on the dwell times and step frequencies of myosin heads under various load 
conditions. We have shown that by a suitable tuning of the internal force and depth of the potential 
well, the theoretically calculated probability p and mean first-exit time /i of the representative Brownian 
particle are in good agreement with their biological counterpart. A second set of experimental data 
concerning the net step number distribution of myosin heads under low load conditions has then been 
exploited to show that the washboard potential used by us is able to reproduce such distribution within 
the mathematical analogy. To the rather small size of the experimental sample should be ascribed the 
3% discrepancy represented by the two net backward displacements predicted by our model. 

Next, the robustness of our model has been tested by inserting 4 different potential profiles in the 
Langcvin equation of motion. The consequently performed calculations have shown that the mean exit- 
time depends on the external force in a way that is essentially insensitive to the chosen potential's profile. 
The chosen profile, instead, has been seen to play an essential role in that it critically relates the depth 
of the potential well to magnitude of the mean exit-time. A finer tuning of the mean exit-time is finally 
achieved by regulating the level of asymmetry of the potential's profile. 
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Tabic 8: The first two columns list the net number of steps performed by myosin heads 
and the corresponding observed frequencies. The remaining four columns show the net step 
distributions, obtained via numerical simulations using the indicated potential profiles all 
possessing Uq = 5 k B T. In all cases F{ = 1.75 pN, F e = 0, xo = L/2, and At = 0.05 ns. 
Other parameter have been chosen as in Figure 2. Variation coefficients and mean-exit times 
Ijl, obtained via Eq. (6), are listed as well. 



\X R /L\ 


Observed frequency 


Potentials' profiles 






Saw-toot 


Parabolic 


Cosine 


Lindner (6 = 2) 


-1 




1.89 


1.92 


1.90 


1.99 





14 


13.61 


13.69 


13.73 


13.81 


1 


21 


20.67 


20.69 


20.59 


20.46 


2 


18 


16.47 


16.42 


16.17 


16.20 


3 


10 


8.63 


8.60 


8.70 


8.64 


4 


3 


3.33 


3.29 


3.38 


3.34 


5 





0.98 


0.95 


1.05 


1.07 


CV 


0.964 


0.970 


0.984 


0.982 


[i (ns) 


1206.048 


1403.504 


2194.007 


2270.677 
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Figure 6: For the saw-toot potential profile with the indicated value of asymmetry 
La, the mean first-exit time fi is plotted as a function of the external force F e . The 
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also indicated on each curve. We have taken Fi = 1.75 pN, while the other parameters 
have been chosen as in Figure 2. 
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